Direct observation of paramagnons in palladium 
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We report an inelastic neutron scattering study of the spin fluctuations in the nearly- ferromagnetic 
element palladium. Dispersive over-damped collective magnetic excitations or "paramagnons" are 
observed up to 128 meV. We analyze our results in terms of a Moriya-Lonzarich-type spin fluctuation 
model and estimate the contribution of the spin fluctuations to the low temperature heat capacity. 
In spite of the paramagnon excitations being relatively strong, their relaxation rates are large. This 
leads to a small contribution to the low-temperature electronic specific heat. 



Nearly ferromagnetic metals are of topical interest be- 
cause their bulk electronic properties can be modified by 
the presence of spin fiuctuations [l|-l5| . Doniach pointed 
out [l| that metals close to ferromagetic order at zero 
temperature should show dispersive overdamped mag- 
netic excitations or "paramagnons" . These should be 
contrasted with the well-defined propagating spin waves 
which occur in an ordered ferromagnetic phase. Over- 
damped modes are still important because they are ex- 
cited with increasing temperature and therefore con- 
tribute to the electronic specific heat [l|-l5|. It has also 
been suggested that they can mediate superconductive 
pairing [a, |7[. In this paper, we report an inelastic neu- 
tron scattering (INS) study of paramagnons in the ele- 
ment palladium. Pd is unique among the paramagnetic 
elements in that it shows a large and temperature depen- 
dent susceptibility Isi. It has one of the highest densities 
of states (DOS) [9l4l2| at the Fermi energy of the d-band 
metals and the measured susceptibility is approximately 
10 times larger than that calculated directly from DOS. 
Thus it is an good system in which to search for para- 
magnons. We find that the paramagnon excitations can 
be observed over a wide range of energy between 25 and 
128 meV in the present experiment. 

Palladium is an face-centered-cubic (FCC) metal with 
lattice parameter a=3.88 A. We studied a 487g single 
crystal of approximately cylindrical shape with a mosaic 
of approximately 1.5 deg full-width-at-half-maximum 
(FWHM). Prior to the experiment the crystal was an- 
nealed at a temperature of 300 C under a vacuum of 
approximately 10~^ torr for 72h to expel hydrogen [l3| . 
Fig. [T] shows the susceptibility of a piece cut from our 
sample compared to a powder standard. Both the refer- 
ence sample and the single crystal used in the experiment 
show an upturn in the susceptibility at low temperatures. 
It is known that even small concentrations of magnetic 
impurities such as Fe can cause such an upturn at low 
temperatures [11] due to paramagnetism of the Fe "giant 
moments" : based on the magnitude of the upturn, we 
estimate the concentration of magnetic impurities to be 
60-90 ppm. 

INS experiments were performed on the MARI instru- 
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FIG. 1: The bulk susceptibility of the Pd single crystal used 
in the present experiment (squares) compared to a standard 
powder. 



ment at the ISIS spallation source. MARI is a low- 
background direct-geometry time-of-flight chopper spec- 
trometer. For the present experiment, we used detec- 
tors located in a single plane henceforth known as the 
scattering plane. The (110) crystal plane was mounted 
coincident with the scattering plane for the present ex- 
periment allowing wavevectors of the type Q = {h, h, £) 
to be investigated. INS probes the E and Q dependence 
of x"(Q,<^). The magnetic cross section is given by 



dV _ 2(7re)^ kf 
dfldE Trg^/ig ki 



F{Q)\ 



x"(Q,M 



1 — exp(— ftw/fcr) 



(1) 



where (77'o)'^=0.2905 barn sr"^, k^ and k/ are the inci- 
dent and final neutron wavevectors and |-F(Q)P is the 
magnetic form factor for a Pd 4d orbital [1^ . Data were 
placed on an absolute scale using a vanadium standard 
and measurements of the acoustic phonons of the sample. 
The relatively large size of the Pd crystal meant incident 
beam was attenuated by absorption and Bragg scattering 
within the sample. We carried out Monte- Carlo simula- 
tions of these effects to account for them in our fitting 
procedure. We use the reciprocal lattice to label wavevec- 
tors Q = ha* + kh* + Ic* . 
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FIG. 2: (Color online) Paramagnon excitations in Pd at 
T—20 K. Data are collected for a single spectrometer setting 
with Ei=71 meV. Excitations are probed over the surface 
of Q — a; space shown in the inset. Phonons are observed 
lower energies hu < 30 meV. The paramagnon scattering can 
be seen near Q=(l,l,-1), which corresponds to 39 meV. The 
black arc is the &j=39 meV contour. The units of the plot 
are mb St~^ meV~^ f.u.~^. 







Pd (T=20K) 




-— 0.04 


- (b) 






^ 


Q=(i,i,i;-i) 


A 




1 0.03 


/ico=38 meV 


- 


m 




f^ f' 




E. 0.02 


- 


A \ 


- 


HI 

S 0.01 

1^" 


■ ^ 

1 i 1 


i 1 i 1 


^y4 





-0.8 






-0.4 




0.0 


0.4 


0.8 




C (r.l. 


u.) 




^ 




1 






'^. 0.03 


- 






- 


=3 


(a) 




1 










1 




> 

^ 0.02 


Q=(-1,-1,C-1) ,, 
/i(B=25 meV 


1 1 1 


1 














w 


' 




"L 




.Q 


1 




4T\ , 




,"-0.01 


-i, 


i"i 


l,# 


1 


* 




.^S 






Li 


iliiT 




1 






-a C. 


Tl 


~v 


-Wfi ' 






"O 


1 


1 I'M 














^ \^ 


T ' 






0.00 





1 


1 



-0.4 -0.2 0.0 0.2 0.4 

C (r.l.u.) 



FIG. 3: Paramagnon excitations at T—20 K. Constant en- 
ergy cuts along (0,0,^) through the (111) reciprocal lattice 
position for (a) hu}=25 meV and (b) fiio=S8 meV. Data were 
collected with Ei=35 and 71 meV respectively. The inte- 
grated proton current delivered to the target during the run 
was 3500 ^Ah. 



In order to avoid any possible complications associated 
with the low-temperature spin freezing |15| , we collected 
data at T = 20 K. Fig. [2] shows data collected with an 
incident energy £'^=71 meV with the [110] direction par- 
allel to ki . The main panel shows the scattering function 
{ki/kf){(Pa/dftdE) plotted as function of wavevector of 
the excitations Q = k^ — ky. Because the data are col- 
lected in a single setting, the energy of the excitations 
probed fvjj = Ei — Ef varies over the figure. The energy 
transfer corresponding to each wavevector can be deter- 
mined from the inset to the figure. At low energies below 
about 30 meV, we observe the highly structured phonon 
scattering. Above the highest phonon energy we observe 
additional scattering near the Q = (1,1,-1) reciprocal 
lattice position. For the present setting this corresponds 
to an energy ficj=39 meV. Fig. EJb) shows a cut directly 
through this position demonstrating that the additional 
scattering is peaked at Q = (1,1,-1). We note that 
there is no observable scattering at the (2,2,-2) recip- 
rocal lattice position with larger |Q|. This is consistent 
with the expected drop in the form factor |i^(Q)| by a 
factor of more than 100 [IJ] expected for magnetic scat- 
tering. Fig. |3l^a) shows that the response near the (111) 
zone center is also present at lower energies. Indeed it is 



peaked at either side of (111). Fig.|4]shows data collected 
at T = 300 K for incident energies E, =35.4, 71, 204.3, 
300 meV. The scattering near the zone center positions of 
(111) and (002) types persists up to the highest energies 
investigated in the present experiment, huj~128 nieV. It 
is interesting to note that width of the response is not 
resolution limited [see Fig. IH^d)] and it broadens with 
increasing energy transfer. 

Magnetic excitations in nearly ferromagnetic and 
weakly magnetic metals are usually interpreted in terms 
of correlated particle- hole pairs (Stoner excitations). Do- 
niach [l| calculated the interacting spin-susceptibility for 
a nearly ferromagnetic metal using the RPA formulation 
of Izuyama et al. |16| , for a single band Hubbard model. 
Moriya [3[ and Lonzarich [J, |5| (ML) developed a gener- 
alized phenomenological form which can be used to de- 
scribe nearly ferromagnetic metals and provides a useful 
input to spin-fluctuation theories. In the ML model the 
imaginary part of generalized susceptibility |l7| takes the 
form. 
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FIG. 4: Paramagnon excitations at T=300 K. Constant en- 
ergy cuts through the (111) and (200) reciprocal lattice points. 
Incident energies used were Ei=35 (a), 71(b), 204(c) , 300 
meV (d). Solid lines are fits to Eq.[T]|4l The additional Une in 
(d) shows the instrumental response to an energy-independent 
sharp response J(Q — r). 
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TABLE I: The results of fit ting the phenomenological ML- 
niodel (Eqs.[2]|4]) to our data [ij. The results [3 of a similar 
analysis for NiaGa are also given. 



where the relaxation rate r{q) is given by 

r{q)^iqX-\q) (3) 

and the wavevcctor-depcndent susceptibility xiQ) — 
x{q, w = 0) is given by: 



X \q)^X ^ + cq^ 



(4) 



The above expansions are valid in the small 9 = |q| limit 
and q is measured from a reciprocal lattice point. Within 
this model the response is characterized by three micro- 
scopic parameters c, 7 and Xj of which only the suscep- 
tibility is temperature dependent. 

We first fitted our r==300 K data to Eqs.[T]|ll The sohd 
lines in Fig.|4]are the results of our fits. Because the bulk 
susceptibility of Pd is well known, we fixed this parameter 
and allowed an overall scale factor to take up any errors 
in the absolute normalization. We found this factor to 
be 0.8 ± 0.2 i.e. unity within the error of the experi- 
ment. The results of fitting the model to the r=300 K 
data are shown in Table ID Within the ML model only 
the susceptibility is expected to vary significantly with 
temperature. The lines in Fig. |3] show the predictions of 
the model using the T = 20 K value of the susceptibility 
from Fig. [TJ It is interesting to note that the response 
at the lowest energy (25 meV) is slightly sharper than 
the model predicts. The doubled peaked structure seen 
in Fig. [SJa) is a consequence of the ML model, because 
r(g) — >■ as g — !■ and does not indicate propagating 
excitations. In spite of being a small q and to model the 
phenomenological ML-model appears to provide a rea- 
sonable global description of the data. The parameters 
in Table. |I] can be estimated from electronic structure 
calculations |10l4l3 |. The calculated values are h'j = 2.1 
ill and 1.1 ^^l meV"^ [H and c= 836 0, 900 [Ilj 
and 925 /U^^ A^ meV 12 1. It is not clear why there is a 
significant discrepancy in the estimation of c. However, 
it has been noted that c is very sensitive to the detailed 
band structure near the Fermi energy |lll . Il2 1 . 

The present results can be compared with those ob- 
tained on NisGa. This material is closer to ferromag- 
netic order at low temperatures and shows a suscepti- 
bility enhancement of about 100 with respect to simple 
band structure calculation [19[. Unfortunately large sin- 
gle crystals of NisGa are not available. Bernhocft et al. 



[l8| carried out an INS study on polycrystallinc mate- 
rial at low energies and found that the response could 
be parameterized in using the model used here. The pa- 
rameters found are shown in Table ID The paramagnon 
excitations in NiaGa are observed at much lower energies 
than in Pd: the maximum energy investigated was a few 
meV. Surprisingly the 7 and c parameters in the two ma- 
terials are the same to within a factor of about 2. The 
energy scale of the spin fluctuations is controlled almost 
entirely by the susceptibility x which is different in the 
two materials. 

As mentioned in the introduction, we would expect 
[2, Isl, la, 120, l2l| that the presence of paramagnons will 
contribute to the low temperature linear specific heat 
C = 'yc'T. Wc may use our phenomenological response to 
estimate the contribution of spin fluctuation to the low- 
temperature specific heat in Pd. Within the ML-model 
y, the electronic specific heat has been estimated to be. 



C = jcT + 6T^ ln(T/T*) 



where 



7c 



A-2 



Anhjc 



ln(f + cxg^), 



(5) 



(6) 



27r4xV(5?i^7), 



and T* ps hy/{kBC'^X'^)- Us- 
c, 7 and X ^-i^d s- cut off wavevec- 



ing our values for 

tor qu at the Brillouin zone boundary ((7bz=1-71 A~^), 
we obtain an estimate of the electronic specific heat of 
7c=5.0±2.7 mJ K^^ mole^^. This should be compared 
with that obtained directly from band structure using 
the the standard relation 7,7 = (7r^/3)fc^A^(e/?) and the 
calculated N{£p)=?>2.7 states atom~^ Ry~^ Q, which 
yields 7c=5.6 mJ K~^ mole~^. We should also note 
that enhancement of the linear specific heat due electron- 
phonon coupling is estimated to lie in the range 28-41% 
22L l23j for Pd and the experimentally determined value 



is 7c=9.42 mJ K""' [2j]. Combining these facts suggests 
that the enhancement in the electronic specific heat due 
to spin fluctuations in Pd is in the range 30-40% which 
is within the uncertainty range of our estimate based on 
the ML-model. Note the estimation based on the ML- 
model is determined only by experimentally measured 
quantities. Thus we have a consistent picture in which 
the observation of strongly spin fluctuations in Pd does 
not lead to a large contribution to the linear specific heat. 
In summary, we have used inelastic neutron scattering 
to measure so-called "paramagnon" excitations in pal- 
ladium. Paramagnons are dispersing overdamped col- 
lective excitations which are present in nearly ferromag- 
netic metals. We observe a dispersing response which 
is strong near the Brillouin zone center and broadens 
in wavevector with increasing energy up to the high- 
est energies investigated, tmj=12S meV. We parameterize 
the observed response and use a Moriya-Lonzarich spin- 
fluctuation model to estimate the low-temperature linear 



specific heat directly from our data. We find that rel- 
atively small enhancement of the specific heat observed 
in Pd is consistent with observed paramagnon spectrum 
which is broad in energy compared to more strongly en- 
hanced systems such as heavy fermions [25|. 
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